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Abstract. Using the multipolar expansion of electrostatic and magnetostatic 
potential energies, we characterize the long-range interactions between two weakly- 
bound diatomic molecules, taking as an example the paramagnetic Er 2 Feshbach 
molecules which were produced recently. Since inside each molecule, individual atoms 
conserve their identity, the intermolecular potential energy can be expanded as the 
sum of pairwise atomic potential energies. In the case of Er 2 Feshbach molecules, 
we show that the interaction between atomic magnetic dipoles gives rise to the usual 
R~^ term of the multipolar expansion, with R the intermolecular distance, but also to 
additional terms scaling as R~^, , and so on. Those terms are due to the interaction 
between effective molecular multipole moments, and are strongly anisotropic with 
respect to the orientation of the molecules. Similarly the atomic pairwise van der 
Waals interaction results in R~^, ... terms in the intermolecular potential 

energy. By calculating the reduced electric-quadrupole moment of erbium ground level 
(J = 6||(52||>I^ = 6) = —1.305 a.u., we also demonstrate that the electric-quadrupole 
interaction energy is negligible with respect to the magnetic-dipole and van der Waals 
interaction energies. The general formalism presented in this article can be applied 
to calculate the long-range potential energy between arbitrary charge distributions 
composed of almost free subsystems. 


1. Introduction 

For a long time few-body effects have been attracting a lot of interest, especially in 
nuclear physics [1], resulting in some striking theoretical predictions like the Efimov 
effect [2]. More recently the tremendous progress for controlling interactions between 
ultracold atoms, has allowed to experimentally confirm those predictions [3]. Indeed the 
signature of triatomic Ehmov states was identihed in an ultracold Bose gas of cesium, 
where the formation of two-body bound states was either forbidden |1] or permitted |5] . 
An essential feature of Efimov bound states is their universality in the following sense 
PmE]: they are characterized by two parameters, the two-body scattering length 
and the three-body parameter [Qj [10], which accounts for all the details of atomic 
interactions. Using ultracold atomic and molecular gases, many extensions of Ehmov’s 
original prediction [2] were then explored or proposed, like four-body [m ESI US] , or 
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more-body bound states [2], the impact of Fermi statistics [T51 [TB] or distinguishable 
particles im, deviation from universality [ISl [181 |T9] , or bound states of dipolar particles 
ED]. Beside Ehmov effect, the theoretical modeling of collisions involving weakly- 
bound dimers revealed the enhanced stability against collisions of molecules composed of 
identical fermions with respect to those composed of bosons [211122]. Besides, few-body 
collisions in the presence of dipole-dipole interactions were also explored [23] . 

In this respect the production of ultracold gases of lanthanide atoms is extremely 
promising [211 ESI EEl [27]. Firstly their strong magnetic dipole moment creates 
anisotropic and long-range dipolar interactions that, unlike electric-dipolar interactions, 
do not need to be induced by an external held. Secondly, erbium and dysprosium 
possess stable bosonic and fermionic isotopes, that were driven to quantum denegeracy 
[2S1 ESI IMl Ell E2]. Despite the absence of hyperhne structure, bosonic erbium and 
dysprosium present dense spectra of Feshbach resonances [33| , which were recently used 
to form very weakly-bound Er 2 molecules, i.e. so-called Feshbach molecules [3l], through 
magneto-association technique. Such Er 2 molecules look like excellent candidates to 
study few-body physics in dipolar systems. 

In many investigations of few-body physics in ultracold gases, atomic interactions 
are described with model potentials, e.g. contact potentials, which is justihed for atoms 
interacting through van der Waals forces. However if the atoms carry a magnetic 
dipole moment, the resulting long-range and anisotropic dipolar interaction requires 
a cautious modeling which takes into accounts the internal structure of the atoms. 
In this article, using the multipolar expansion in inverse powers of the intermolecular 
distance R [35l [36], we characterize the long-range interactions between two weakly- 
bound diatomic molecules. We focus on the regime where the two molecules are 
approaching each other, that is to say when the intermolecular distance is larger than 
the mean interatomic distance inside each molecule. Assuming that individual atoms 
keep their own identity within each molecule, we expand the intermolecular potential 
energy as the sum of pairwise atomic interaction energies. Taking the example of two 
Er 2 Feshbach molecules, we show that, when expressed in the coordinate system of 
the molecule-molecule complex, the total dipole-dipole and van der Waals interactions 
between all atom pairs are both expressed as a sum of terms proportional to inverse 
powers of R, and involving effective molecular multipole moments. These terms, which 
are absent in the usual multipolar expansion, are strongly anisotropic with respect to 
the orientation of the two molecules. In addition, by calculating the electric-quadrupole 
moment of erbium ground level we show that the total quadrupolar interaction, which 
would in principle give rise to another series of terms, is actually much smaller than 
the dipolar and van der Waals interactions. We will calculate adiabatic potential-energy 
curves that could be used in a future work to study Er 2 -Er 2 collisions in the ultracold 
regime. 

The article is outlined as follows. In Section E] we present the general formalism 
giving the hrst-order and second-order energy corrections between arbitrary weakly- 
bound charge distributions. This formalism appears as a generalization of the usual 


Four-body long-range interactions between ultracold weakly-bound diatomic molecules 3 



Figure 1. (Color online) (a) Vectors, (b) distances and angles describing the position 
of the weakly-bound diatomic molecules A (composed of atoms 1 and 2) and B 
(composed of atoms 3 and 4) in the space-fixed frame XYZ with quantization axis Z. 
The relevant coordinates describing a pair of atoms from A and B is drawn, with the 
example of atoms 1 and 3. For sake of clarity identical atoms are considered, and the 
azimuthal angles are not represented. 


multipolar expansion. Then in Section |3] we consider the example of Er 2 Feshbach 
molecules, focusing on their magnetic-dipole and van der Waals interactions, and 
discussing also their electric-quadrupole interactions. With arguments based on the 
characteristic lengths associated with the multipolar expansion, we show that the 
anisotropic terms due to effective molecular multipole moments, are likely to play 
an important role in the Er 2 -Er 2 collisional dynamics at ultralow energies. Section 
m contains concluding remarks. 

2. Potential energy between distant weakly-bound molecules 

We consider two molecules denoted A and B. Molecule A is composed of atoms 1 and 2, 
and molecule B of atoms 3 and 4. The orientation of the interatomic axes of A and B in 
the space-fixed (SF) frame XYZ, Z being the quantization axis, are characterized by 
the vectors Ra = Ri — R2 and R^ = R3 — R4, with spherical coordinates {Ra, ©a, '^’a) 
and {Rb,Qb,^b) respectively. The orientation of the intermolecular axis, which joins 
the centers of mass of A and B is given by the vector R = [R, 0, $) (see Fig. [1]). 

2.1. First-order correction 

When the distance R goes to infinity, molecules A and B are independent from each 
other. Their quantum states, called |A) and \B), are characterized by the zeroth- 
order energies Ea and Eb respectively. As R decreases, A and B start to interact 
through electrostatic and/or magnetostatic forces. We assume that within the weakly- 
bound molecules, each individual atom keeps its identity - namely that the exchange 
is neglected - so that atoms interact through electrostatic and/or magnetostatic forces 
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Table 1. Coordinates of the vectors Ri, R 2 , R 3 and R 4 , as well as the geometric 
factors gi, 772 , Vs and 774 appearing in the intermolecular potential energy (HU, given 
as functions of the coordinates of R ,4 and Rb, and of the atomic masses M. 2 : 
Ms and M 4 . 


atom 

1 

2 

3 

4 




Ri 

R 2 

R 3 

R 4 


M 2 

M 1 +M 2 

■Ml 

M1+M2 

Mi 

.M3+.M4 

Ms 

Ms-hMi 


Oj 

Ra 01 = 0 a 
Ra 02 = TT - 0 a 
Rb 03 = 0 S 
Rb 04 = TT — ©B 


4>i, _ 

<l>i = 

4*2 = 4“ 

^3 = 

$4 = $5 + TT 


Vi, Vj 
Vi = k 
V2 = - 

V3 = k 

Va = - 


M2 

M1+M2 

Ml 

.M1+.M2 

Mi 

M 3 +M 4 

Ms 

Ms-hMi 


as well. The total potential energy of the complex V^ot can therefore be expanded as 
the sum of pairwise atomic energies, Vtot = V12 + V34 + SLi Si=3 work, 

assuming that the hrst two terms V 12 and V34 are part of the unperturbed energies Ea 
and Eb, we focus on the intermolecular potential energy 

y(RA,RB,R)=EE^o(R7i)> (1) 

i=lj =3 

where Rjj is the vector pointing from atoms i E A to j ^ B. 

In Eq. ([1]) the interatomic energy Vij is given by the usual multipolar expansion in 
the SF frame 

+00 

r«(R.,)=E E ( 2 ) 

^ij — 0 TTi-ij — 


where the factor -Fq only depends on the atomic multipole moments. 


FeijrriB - Fo E (“1)^" 

q,q=o 

+ii +b 


24 

2ii 


1/2 


X E E 

mi=—£i mj=—ij 


(3) 


with Eq = l/47reo (po/47r) for electrostatic (magnetostatic) interactions, cq and po the 
permitivity and permeability of the vacuum, (:) a binomial coefficient, = 

{iimi£jmj\£i£jiijmij) a Clebsch-Gordan coefficient, and Qqmi {Qe^mj) the multipole 
moment of atom i (j), expressed in a coordinate system (CS) of the SF frame and 
centered on atom i [j). The multipolar expansion can be applied if all interatomic 
distances are larger than the so-called LeRoy radius 1311 . SO that their electronic clouds 
do not overlap. 

In Eq. ([ 2 ]), G(,^^rnij is a purely geometric factor involving the Racah spherical 
harmonics (0^^, 4)*^). 
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It depends on the relative coordinates of atoms i and j, which are not compatible to 
the CS dehned in Figured! In order to express in this CS we use the following 

relation 


Rjj = R - (Rj - Kj) = R - IJji, 


(5) 

where Rj (Rj) are the vectors joining the center of mass of molecule A (B) to the atom 
i (j). Those vectors are colinear to Ra (Rb) and their coordinates are given in Table 
dl The coordinates of U^j are called {Uji, Eji, 'hjj). 

To transform Eq. (jT]), we expand the spherical harmonics for the vector r = r' —r" = 
as function of those for the vectors r' = R and r" = U^j. Setting r = (r, 9, 0) (and 
similarly for r' and r"), we apply (see Ref. [38!, C^h. 5, Eq. (36), p. 167 T 


Ck,{0,<P) ^ 


{2k' + l)\ 


1/2 




-.l+Zc 


{2k'')\{2k + 1)\ I {F)k'+i 


k',k"=0 

+/c^ -\-k" 

X i; i; c^w'twCi.'v'(o".<t'")Ckv(o',4‘') 


( 6 ) 


q'=-k' q" = -k" 

which is valid for r" < r'. We obtain for Eq. (0!) 


+ 00 


(R-p, 


= E (-1)^ 


(2A + 1)! 


1/2 TT^ (-ij 

1 ji 


\=L. 


X2A-2%)!(24 + l)! 




+A 


X 


/j.=— A 








which is then valid for Uji < R. Since 0 < Uji < Ri + Rj, Eq. ([7]) is applicable provided 
that 


Ri-\- Ri < R. 


( 8 ) 


We will come back to this criterion later on. 

Equation d?]) represents the hrst main result of our approach: the pairwise atomic 
potential energy has been transformed from a sum of terms in the relative CS of the 
interacting atoms (see Eq. (jTj)), into a sum of terms proportional to inverse powers of 
the intermolecular distance R (see Eq. ([7])). The price to pay is the emergence, for 
each couple of atomic multipole moments of an inhnite sum of terms scaling as 

^-i-A term is anisotropic due to the Racah spherical harmonics C'A^(0,‘h). To 

express the energy as a function of Rj and Rj, we apply a second transformation (see 


Ref. 


, Ch. 5, Eq. (35)) to r = Ujj, r' = Rj and r" = Rj 


'r’^Ckq{6,(j)) — ^ dk'+k'',k{—^Y 

k',k"=0 

+k' +k" 


{2k)\ 


1/2 


{2k')\{2k'‘ 


(r'Yir'T 


q' = — k' q" = — k" 


(9) 


f In [38] the relation is given in terms of normalized spherical harmonics Ykq] we transform it using 
Ykq{e, 4>) = v'(2fc + l)/ 47 r X Ckq{e, b) 
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which is valid for all r' and r". This gives the final expression for 


G 


+ 00 ^ — ^ij 

A —£ij A^jAj—0 


- 1 ) 


A+Aj 


2A + 1 

2iij + 1 


2 Aj T 2 A} 

2A,; 


1/2 


X 42^ 1: 


jl=—\ 

+Ai +Aj 


X 


— Aj flj — \j 


( 10 ) 


In addition to the snm over A, we get two snms over Aj and \j which, dne to their 
dependence (and similarly for j), can be viewed as effective mnltipole 
moments describing the position of atoms i and j in the CS associated with A and B 
respectively. It is worthwile mentioning that for i?,, Rj R, namely when the size of 
the molecnles is negligible with respect to their mntnal distance, the snms rednce to 
Aj = Aj = 0 and to A = ^ij, and so we recover the nsnal mnltipolar expansion. 

Finally, when adding np all pairwise atomic contribntions, we can replace Aj by 
Aa and Aj and As, as the position of atoms 1 and 3 in the CS of A and B gives the 
orientation of the interatomic axes of A and B respectively. For atoms 2 and 4, the 
orientation is opposite to that of the corresponding interatomic axes, and so we nse 
Cfcq( 7 r — 6 ^, 0 + vr) = {—lYCkq{0^ 0). Moreover replacing Gj by ^ for simplicity, and nsing 
Ckg{0,(l)) = {—iyCk-q{0,(l)) and the particnlar expression of Ckfqik"q" for k = k' ± k" 
|38] . we obtain for the intermolecnlar potential energy of Eq. o 


+00 +00 \—l 

E(Ra,Rs,R)=Eo^ 5 : ^ 

d\A+\B+(;>^ (“1) 

£=o x=e Xa,\b=o 

+A 


E)Aa td^b 


X 1 : ^(A + f,)!(A-f,)!xC;^(e,4) 

/i=—A 


X 


£IA=-^A Ais=-As \/+ Fa)K^A — Fa)K^B + Fb)K^B — FbV- 


i=i j=3 q,q=0 

X E E 


M A+Ms +Ej rrij 


\/(£i + - m,)! 


( 11 ) 


where the geometric factors rji and rjj are given in Table [H The Kronecker symbol in 
the last line of Eq. ffTTj) is obtained by combining the Clebsch-Gordan coefficients of 
Eqs. (El) and (ITOll . Since the condition fiA + I^b + fni -\- mj = ja mnst be satisfied for all 
pairs (i, j), it imposes that mi = m 2 and m 3 = mi. 

As already mentioned, Eq. oa can be viewed as the interaction between effective 
mnltipole moments of ranks Aa and As, describing the orientation of the interatomic axes 
of molecnles A and B. Eqnation (ITT|l comes ont as the nsnal snm of terms proportional 
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to inverse powers of the intermolecular distance and to angular factors <h), 

which account for the anisotropy of long-range interactions. Since —1 < rji^rij < 1, 
Eq. (IH]) can be rewritten as a condition of validity for Eq. (mi 

max(|)),|, \rh\)RA + max(|t; 3 |, lihD^B < R- (12) 

In the honionuclear case, e.g. Adi = A42, ma.x{\T]i\,\'r] 2 \) = 1/2. On the contrary if 
Ail > A42 or Ail -C A^ 2 , then max(|? 7 i|, \r] 2 \) = 1. 


2.2. Second-order correction 


Calculating the second-order correction to the intermolecular potential-energy ([T]) is 
equivalent to calculating the matrix elements between unperturbed states of the second- 
order operator 


IE(R^,Rb,R) 


E 

{A" 


V\A"B''){A"B''\V 
Ea" — Ea + Eb" — Eb 


(13) 


where \A'') and \B'') formally denote the states of molecules A and B which are 
coupled to I A) and \B) by V. The sum is performed for all possible pairs of states 
{\A”),\B")) excluding the case where both \A") = jA) and \B") = \B). Because the 
unperturbed states |A) and \B) correspond to weakly-bound molecules, we can suppose 
that the states \A"B") such that {AB\V\A"B") ^ 0 also correspond to weakly-bound 
molecules, but near different atomic dissociation limits. Therefore in Eq. (1T3|1 the energy 
differences Ea" — Ea and Eb" — Eb can be replaced by the energy differences between 
the corresponding atomic dissociation limits. We can replace the sum over the molecular 
states \A'') and \B'') by a sum over states of the separated atoms |1"), |2"), |3") and 
|4"). This assumption implies that the geometric factors given by Eq. ffTOj) will 

be taken out of the sum over atomic states, which gives for Eq. (IT^ 


+00 






w(Rx,Rb,r)= - E E E E E 

j,j'=3 liAl^iji=Q mij=-lij 

J2'/3", 4'')(1'', 2", 3", 4"|F« 


X E 


Ai" -|- A2" -l- A3// -|- A4 


(14) 


where A^// = E^// — E^ is the excitation energy of atom A; (A; = 1 to 4). The sum over 
atomic states excludes the case [A:") = \k) for all the atoms at the same time. 

Applying Eq. ([Tl]) with the particular form of and given by Eqs. ([3]) 

and flTUD would be inconvenient, as it would not yield terms with irreducible tensors, 
and so would prevent from using the Wigner-Eckart theorem for the matrix elements of 
W. Insted we introduce coupled multipole moments of rank kA associated with atoms 
i and i' (see e.g. iMBniiiii) 


mi=-li m ■/=-£■/ 


( 15 ) 
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and the same for j, j' and B. Doing similar transformations for Racah spherical 
harmonics (see Appendix A for details), we obtain the hnal expression for W 


ir(r^,Rb,r)= -f^ y. E 


,^A^+A'. \g+\' 


,6 


TD''A I ^'A p 
Ih A Ih 


B 


eXXA^B ^'A'A' A' 


'^A“l“'^s“l"^^d'^A/i+As+£,A ^2+A+A' 


X 


2A + 1 \ / 2Aa + 2Ab \ / 2A' + 1 \ / 2A'. + 2A', 


-,1/2 


B 


2i+l 


2A, 


2f + 1 


2Ah 


X y [(A^ + Ab)(A(4 + 

kK.TTK.Af^'B 



' Aa 

As 

Aa + As 


Aa + As 

A 

^ ] 


A'a 

As 

A'a + A's 


a: 4 + A's 

A' 

f 


Ka 

Kb 

71 

> 


- 

K 

k J 




X 


E c 

qQCTQA Bb 


Ayi 0 A'^ o'-'Afl 0 A'g o'" AO A' 0 

CS„C-(9, tjC* „ (0^, 4.4)C- „ (Ob, 4b) 


f^AQAl^B Qb 


iiiil 
1/2 




X 


2i 

2L 


2 f 

^21, 


ii ij 


X Yc^^ 


E 


m Y [kAkB^^^ { ty I' 

kA ks k 


kAkg 

Q" 


{£i,£ii)kAqA ~^{^j,£ji)kBqB 


449B A” + ^2" + A3B + A4/, ’ 


(16) 


This eqnation contains several snms: over all the qnantnm states |1"), |2"), |3") and 
|4") of the fonr atoms, on the atoms themselves {i,i' = 1,2 and j,j' = 3,4), and over 
tensor-operator ranks and components. In this respect the Latin letters correspond to 
the atomic mnltipole moments, and the Greek ones to the effective molecnlar mnltipole 
moments. The letters i and A characterize nnconpled tensor operators, whereas k, k 
and 71 characterize conpled ones (see below), and q, g and a their components. The 
nnprimed (primed) nnconpled tensor ranks come from the hrst (second) call of the 
mnltipolar operator V in the second-order correction (see Eq. flT^ h 

Eqnation (ITHll shows that the second-order mnltipolar interaction is based on 
bnilding blocks which are tensor operators, i.e. the atomic mnltipole moments ii, ij, 
iii and ij>, and the effective molecnlar ones Aa, As, Aa' and A^/. The /^-dependence 
of the operator W is obtained by adding those bnilding blocks, namely . The 

angnlar dependence of W is associated with conpled tensors whose ranks are obtained 
by adding np the bnilding blocks in the sense of angnlar momentnm theory (see Table 

ED. 


3. Example: interactions between Er 2 Feshbach molecules 

In a recent experiment [3l], an nltracold gas of bosonic ^®®Er atoms (with vanishing 
nnclear spin / = 0) was prodnced in the lowest Zeeman snblevel \J = Q,Mj = —6) of 
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Table 2. Mathematical and physical definitions of the ranks of the coupled tensors 
appearing in Eq. (HU). They are constructed by vector addition in the sense of angular- 
momentum theory. 


tensor rank 

physical quantity 

kjA 

kB — + ^j' 

l^A = + ^A' 

Kb = + As' 

K = X + X’ 

coupled multipole moment of atoms {i,i') 
coupled multipole moment of atoms {j, j') 
coupled effectivemultipole moment of molecule A 
coupled effectivemultipole moment of molecule B 
coupled tensor for the orientation of the intermolecular axis 

71" — (Aa + A^') -|- [Xb + Xb') — ka + Kb 
k = Tr-\-K = i-\-i' 


the atomic ground level [Xe]4/^^6s^ ^Hgi / = 0. A magnetic-field ramp was applied in 
order to transfer pairs of free atoms into a weakly-bound molecular level, thus creating 
a so-called Feshbach molecule. For molecule A, such a level called |n^) can be expanded 
in the general form 

NaAIn^ 

where Mj^ and Mj^ are the projections of total electronic angular momentum of 
resp. atoms 1 and 2 on the magnetic-field axis Z, Na and are the angular 

momentum and its projection associated with the rotation of the interatomic axis of 
molecule A, Ra is the distance between atoms 1 and 2, and 

multi-channel radial wave function describing the rovibrational motion of the molecule. 
The couplings between the different channels \Mj^Mj^NaM]si^) are due to the magnetic- 
dipole and van der Waals interactions between two erbium atoms. Since the entrance 
open channel corresponds to the atoms in the lowest Zeeman sublevel colliding in s 
wave, i.e. Mj^ = Mj^ = —J = —6, Na = = 0, the allowed channels in Eq. flTTl) 

are such that [3l]: (i) Na is even and (ii) Mj^ Mj^ = —2J = —12. The same 

selection rules apply for molecule B. 

3.1. Magnetic-dipole interaction 

Each erbium atom carries a permanent magnetic dipole moment equal to —fiBdjJ, J 
being the electronic angular momentum (J = 6), with /i^ the Bohr magneton and 
gj = 1.16683 ~ 7/6 the Lande ^f-factor of erbium ground level ^Hq. Following Eq. (1TT|1 
the first-order interaction is such that £* = = 1, £ = 2 and Fq = /xo/dTr. Due to 

that interaction two Er2 Feshbach molecules in levels \va) and \vb) colliding in the 
partial wave L and Z-projection can undergo elastic or inelastic scattering towards 
\v'j(v'^L'M'^) . The matrix element of the magnetic-dipole interaction V^d is then 

{vAVBL'M'^lVj^dlvAVBLMi) 
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’ .. .'i 


X E E 
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^(1 + M',+M,,)! 
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X E E 
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^{1 + MX - Mj,)\{l - M'j^ + Mj,)\^ 

r-\-oo f 

E E L dR,RY Xm'j M'j (-^ a ) XMj^Mj^NaMn^ (Ra) 




X E E 


r+oo 


dRs R^B Xm'j M'j N'gM'j, (Rb) XMj^Mj^NBM^g (Rb) 


X ^(A + Ml- Mi)\{X -Ml + M')!^ 


2 -^' + 1 ^Lo r-LML 


2L + 1 


^LO n 

A,Ml 


-^'a^'na 


X 


X 


I - n‘ 

I2Na + 1 '-'NaO\a(FNaMma>^a,M'nj^-Mma 

'iN'A + 1 + M„J 


2 Nb + 1 


]\fi 0 N'„M'„ 

n B Ng 

'^NBO\BO'^NBMNg\B,M'ffg -Mnb 


+ 1 \/(Ab + - M„J!(Ab - 


( 18 ) 


Contrary to Eq. flTT|) the terms arising from the pnrely atomic part of the interaction are 
written in the first five lines of Eq. (1T8|1 . The atomic dipole moment is expressed using 

-4 J' M' 

the Wigner-Eckart theorem {J'MjJQimi= -iaBgjCjj^XXmAd'\\J\\J)/x'‘^J + 1 
with {J'\\J\\J) = J{J + 1)(2J + 1) the reduced electronic angular momentum. For 
homonuclear molecules, only even values of Aa and \b are possible, since gi = —772 = 
7^3 = —7^4 = 1/2. All the tensor components of Eq. (fTTj) are replaced by their only 
possible value, i.e. mi = M'j^ — Mj^, and similarly for atoms 2, 3 and 4, /i = Ml — M/, 
fiA = ^Na ~^Nai similarly for B. The relationship between the tensor components 
established in Eq. flTT]) imposes the following selections rules 


M'l + dR^A 4“ ^'nb R R dTjg = Ml + M^Ia + Fl^g + Mj^ + Mj^ (19) 

M'l + AfE + ^'nb R A^J2 R MX = Ml + Mj^a + M^g + Mj^ + Mj^ (20) 

M'l + Af^ + ^'nb a ^'ji + ATj^ = Ml + M^Ia + Afjv^ + Mj^ + Mj^ (21) 
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M'l + + M'nb + ^j2 + ~ + ^Na + ^Nb + ^J2 + ( 22 ) 

which implies 


M',, - M'j^ = Af,, - Mj, (23) 

A4 - M'j^ = Af* - A/j.. (24) 

Finally in level I^a) (see Eq. (1T7|1 ). the rotation of the interatomic and intermolecular 
axes is described by spherical harmonics (©a,' f’yi), YnbMnb^^b^^b) and 

"h), and the integral of products of three spherical harmonics is used to obtain 
the last three lines of Eq. ffTSj) . 

Equation ffTSjl consists in a series of inverse powers of the intermolecular distance 
R. The leading term, which scales as appears for Aa = As = 0 and A = 2. It 
couples the channels characterized by the same Nb and M)v^, but by possibly 

different Ml and MJ^.. Taking M'^ = Ml and Mj.. = Mj^. yields the usual two-body 
dipole-dipole interaction (see Eq. (|2])) between magnetic moments d^A and dyg such that 


dvA--FBgj + ^- 72 ) 

with 

f+°° / \ 2 


(25) 


(26) 


and similarly for B, 3 and 4. In Ref. [M] the Er 2 -Er 2 magnetic-dipole interaction 
energy was calculated using such a two-body expression with experimental values of 
dyj^ and . Indeed evaluating quantitatively each term of Eq. flTSl) requires to know 
in details the nature of the Feshbach states, namely the functions (-^a) 

and (-^b)) which is not possible for the moment in Er 2 . 

As a consequence we consider a simple model, where the two molecules are in the 
same state va = vb = v imposing a single-channel condition. As each molecule is made 
of two atoms in the lowest Zeeman sublevel and colliding in s wave, we assume that the 
resulting molecules are in the d-wave bound level, Mj^. = —J = —6, Na = Nb = N = 2 
and = 0, so that the R~^... terms appear in Eq. (IT^ . The mean 

interatomic distance inside each molecule is (uaI-RaI'^a) = {vb\Rb\vb) = Rq. We also 
assume that Eq. flTSl) does not couple different molecular states, i.e. v'a = va and 
v'b = Vb, but that it couples different partial waves L and L'. The molecules are 
assumed to collide in the s wave, which is the case in the temperature range of Ref. |34j . 

Combined with the selection rules (1T^ - (|2T|) applied for v'a = va and v'b = vb, the 
s-wave condition imposes that M'j^ = Ml = 0 for all states. Moreover the Clebsch- 
Gordan coefficients of Eq. flTSl) impose Aa = As = 0, 2 and 2N = 4, and so A = 2, 
4, ..., 4A^ -|- 2 = 10. To evaluate the importance of each term, we compute adiabatic 
potential-energy curves, obtained after diagonalization of the Hamitonian 


Hi{R) 


n^b 


+ Rmd(R) 


(27) 
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where A^red is a Er2-Er2 reduced mass, and L the dimensionless angular momentum 
associated with the collision, and l^d is the magnetic-dipole term given by Eq. flTSl) , in 
the basis spanned by L (note that = 0), for various intermolecular distances R. 

The results of such calculations are presented on Fig. [2], when the lowest adiabatic 
potential-energy curve is plotted in different situations. To highlight the influence of 
the different terms, we truncate the sum on A in Eq. ffTSj) up to Amax, which ranges 
from 2, corresponding to the R~^ term, to its largest possible value 10, corresponding 
to the R~^^ term. Partial waves from L = 0 to Tmax = 30 are included in the basis 
for a proper convergence. We assume that (R^) ~ (Ra)^^ = Rq^ (and similarly for 
B), which can be justified by the expected strong localization of the vibrational wave 
function around the outer classical turning point of the corresponding potential-energy 
curve. We take Rq = 80 bohrs as a typical value for d-wave resonances observed in |34j . 
Strictly speaking our calculation is not valid for R < Rq and so the lower bound of the 
X axis in Fig. [2] should be i?o- But we extend it to a slightly shorter value since Rq is 
only an estimate. 

Figure [2] shows that higher-order R~"' effective-multipole terms {n > 3) get more 
and more important as R decreases. At R = Rq all the terms bring similar contributions 
to the potential energy. More surprisingly, even a.t R = 120 bohrs = 1.5 x fhe R~^ 
term (u max = 3) only accounts for two thirds of the total energy. Because the range of 
energy on Fig. [2], of a few mK, that is tens of MHz, is larger than the typical energy 
spacing between neighboring Feshbach levels, couplings with those levels should be taken 
into account in order to give more accurate predictions. 

After the magnetic-dipole interaction, the next term of the multipolar expansion 
between two erbium atoms is the electric-quadrupole interaction. We have calculated 
the reduced quadrupole moment of erbium ground level (J = 6||Q2||d = 6) using a 
Dirac-Hartree-Fock method, and found -1.305 a.u. We have evaluated its impact on the 
intermolecular interaction by diagonalizing the Hamiltonian 

+ ( 28 ) 

where the electric-quadrupole energy Wq is obtained by setting ii = ij = 2, i = 4 
and Fq = l/dvreo in Eq. flTTl) . The latter shows that the electric-quadrupole interaction 
consists in repulsive contributions scaling from R~^ to R~^^. Figure [2] shows that the 
influence of the quadrupole interaction is visible for i? < 90 bohr, a region where the 
van der Waals interaction will be actually dominant 

3.2. Van der Waals interaction 

The next term of the pairwise atomic multipolar expansion comes from the van der 
Waals (vdW) interaction, proportional to Ri^^■ In Ref. |13] we have shown that 
the Er-Er vdW interaction is mostly isotropic, and characterized by a coefficient 
= 1760 a.u.. Therefore in this section we calculate the second-order electric- 
dipole interaction between two weakly-bound diatomic molecules whose atoms interact 
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Intermolecular distance R (bohr) 

Figure 2. (Color online) Lowest adiabatic potential-energy curves characterizing the 
magnetic-dipole interaction between two ^®®Er 2 Feshbach molecules in the same level 
and colliding in the s wave [see Eqs. (HID and (EZD]. in Eq. (HU the sum on inverse 
powers n = A -I- 1 of the intermolecular distance R is either stopped at Umax = 3 (solid 
line), n„iax = 5 (dotted line), rimax = 7 (dashed line), rimax = 9 (dash-dotted line), or 
the largest possible value rimax = 11 (dash-dot-dotted line). The curves with squares 
also accounts for the electric-quadrupole interaction (see Eq. (1281) 1. The main panel 
is in log-log scale, and the inset is a linear-scale zoom on the small-distance region. 
Other parameters are: Rq = 80 bohr and Lmax = 30. 


through an isotropic vdW term. Such a calculation is applicable for Er2-Er2 interactions, 
but also for molecules made of alkali-metal or alkaline-earth-metal atoms, for which vdW 
is indeed the strongest interaction. 

The vdW energy is a second-order correction due to the electric-dipole interaction, 
= ij = ii' = £j' = I, i = F = 2 and Fq = l/dvrco in Eq. flTB]) . while the isotropy 
results from = ks = k = 0. Pointing out that k = 0 implies n = k we obtain for 
states va = vb = V 


{vvL't) I W^Aw\vvLt)) 

+ CXD 

E 


^Er—Er 
'-' 6,000 


6(47r£o)^-R® 


A A' 

E E ^ A ^-t-A -t- 2 , A ^ A +A -t- 2 , A^ 

A,A '=2 Aa,As= 0 , 2 ,... A'^,A) 3 = 0 , 2 ,... 


'W) 




X 




















Four-body long-range interactions between ultracold weakly-bound diatomic moleculeslA 


X 


E 


Aa 

AA 


Ka 


As Aa + As 

Ai, AUA's 

Kb 


K 


Aa + As A^ + X'b k 
A' A 2 


pK,B^ pK,0 pnO 

^ ‘-'AaOA^o'-'AsOA'^o'-'AOA'O'-'kaOkbOI 


12L' + 1 
2L + 1 




LO \ 

L'OkO^A^O«:aO^A^OkbO i 


where we used 




Er—Er _ 

6,000 




6,000 


= - E 
2 ^ 
" -•// -•/> 


O" O" 

^(£i=l,£i=l)00^(£,=l,£,= 


1)00 


Aj// + Aj 


(30) 


and 
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with {:::} a Wigner 6 -j symbol, and where (i?^ ^ '^’®) is the interatomic distance at the 

corresponding power in molecules A and B, averaged over the vibrational wave function 
of state |n). Equation fl2^ is a series of terms proportional to R~^. The leading term, 
which scales as R~^, comes out when Aa = As = A (4 = A'^ = 0 , and so ka = kb = k = 0. 
It is thus a fully isotropic term (L = L') equal to —ACf^^QQ^/R^. The next terms scale as 
R~^, R~^^, .... Unlike the first-order expression flTSD . n goes a priori to infinity, since it is 
not limited by the angular selection rules. The bounds in the sums over ka, kb and k are 
not explicitly specihed, as they come from several conditions. The 9-j symbol of Eq. 


imposes |Aa — A( 4 | < ka < Xa + X'^, whereas the Clebsch-Gordan coefficient of the last 
line imposes 0 < ka < 2A and ka even. The most restrictive conditions will indeed 
apply, namely max(|AA — A(4|,0) < ka < min(AA + X'j^,2N), ka even. The conditions 
are similar for kb, while for k we have: |Aa + A^ — A (4 — A'^| < k < (Aa + A^ + A (4 + A^), 
I^A — i^b\ < < («^A + i^b), \L — L'\ < k < {L -\- L'), and k even. 

The convergence of the R~^ series in Eq. fl2^ is addressed on figure [3], where 
we plot the lowest adiabatic PEC obtained after diagonalization of the hamiltonian 
H 2 {R) = H[{R) + lUvdw(-R) [see Eqs. fl28|) and fl29l) ] including the magnetic-dipole, 
electric-quadrupole and van der Waals interactions, and where again we assume 
^ j^^a+>^a+^b+>^b^ Contrary to the hrst-order correction, convergence 
will require higher terms as the ratio Rq/R decrease. At i? = 100 bohr, that is to say 
Rq/R = 4/5, convergence is reached for Umax = 40, whereas for R < Rq, the condition 
r" < r' to apply Eq. ([ 6 ]) indicates that convergence cannot be achieved. In addition, we 
see that the van der Waals energy is significantly larger than the first-order energies on 
the left part of Fig. [3l at i? = 120 bohr it represents 73 % of the total potential energy. 


3.3. Ultracold collisions and characteristic lengths 

In order to estimate the role played in collisions at ultralow energies by the 
additional terms of the multipolar expansion, it is instructive to calculate the so-called 
characteristic length associated with each of those terms |4n ITS] . In our case, a given 
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Intermolecular distance R (bohr) 

Figure 3. (Color online) Lowest adiabatic potential-energy curves characterizing 
the magnetic-dipole, electric-quadrupole and van der Waals interactions between two 
^®®Er 2 Feshbach molecules in the same level and colliding in s wave. In Eq. (|2^ 
the sum on Xa, Xb, X'a and X'g is stopped at different values of Umax, where 
n = 6 -|-A/i-|-Ab-|-A( 4 -|-A^, namely rimax = 6 (solid line), Umax = 10 (dotted 
line), Umax = 20 (dashed line), Amax = 40 (dash-dotted line) and Umax = 50 (open 
circles). For comparison the lowest adiabatic PEC only accounting for the first-order 
magnetic-dipole and electric-quadrupole interactions is shown again (dash-dot-dotted 
line). Other parameters: i?o = 80 bohr and Lmax = 30. 


term will cause significant reflection of the incoming scattering wave function, if its 
characteristic length is larger than the typical size Rq of each molecule. 

We start with discussing the influence of the magnetic-dipole interaction. We 
assume that each term of Eq. fflSj) is described by a coefficient 


{2yiBgjjf , (32) 

where we take the largest possible magnetic moments = d^g = —2yt,BgjJ for 
molecules A and B. The characteristic length associated with the term Cn/Rl^ is [^H5] 


f 2MredCn \ 

2 \ ) 


(33) 


By inserting fl52]) into fl5^ . we can express all the characteristic lengths as functions of 
Os = Aij-edCs/h"^ 


dn 

Rq 


Ro 


1 

n —2 


(34) 


which gives (a^/Ro) = {as/RoY^^, {a^/Ro) = {as/RoY^^, and so on. 

Equation fl3Tl) shows that if > Ro, then as > > aj > ... > Rq, for all n, and 

vice versa. So if the R~^ term signihcantly influences the ultracold dynamics, the R~^, 
R~^, ... terms associated with effective molecular multipole moments of higher ranks 
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scaled molecular size Rq / Sp 


Figure 4 . (Color online) Characteristic length a„ as function of the size Rq of the 
molecules. Both quantities are rescaled with respect to their “parent” characteristic 
length, namely ap =3 for magnetic-dipole and ap=Q for van der Waals interaction (see 

Eq. (1551)1. 


are also likely to do so. This is indeed the case in our present Er 2 -Er 2 study where 
Rq = 80 bohr, 03 = 790 bohr, 05 = 172 bohr, 07 = 126 bohr, etc. This reasoning 
can be generalized to all pairwise atomic interactions. For a “parent” term scaling as 
CpjBP and associated with the characteristic length Op (see Eq. ([33])), the related term 
Cn/R^ = Rq~^Cp/K^ is associated with the characteristic length 


Rq 


Rc 


P-2 
n —2 


(35) 


for n > p. Considering the van der Waals interaction p = 6 with Cq = A x = 

7040 a.u., we obtain oq = 128 a.u., Og = 109 a.u., Oio = 101 a.u., etc. The terms coming 
from the van der Waals interactions are also likely to play a crucial role. 

Finally we investigate the influence of the size of individual molecules, which, for 
halo molecular states close to a Feshbach resonance, is equivalent to the atom-atom 
scattering length |16] . To that end we put the Rq term in the right-hand side of Eq. 
and rescale both sides with respect to Op 



Figure m shows this quantity as function of the scaled molecular size, for the three terms 
? 7 ,=p-|- 2 , p-|-4 and p-\-6 associated with the magnetic-dipole p = 3 and van der Waals 
interactions p = 6. When the size of the molecule is much smaller than the parent 
characteristic length {Rq <C Op), the scaled characteristic lengths a„/ap vanish and the 
important interactions are the dipole-dipole interaction p = 3 and the van der Waals 
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interaction p = 6. But when the molecular size increases up to the parent characteristic 
length (i?o ~ Op) 5 On —t Op and all the higher molecular effective multipole moments 
become signihcant in the overall interaction and hence in the ultracold dynamics. 

4. Concluding remarks 

In this article, we present the general formalism to characterize the long-range 
interactions between two arbitrary charge distributions, each composed of almost free 
subsystems, which we apply to the case of weakly-bound diatomic molecules. To that 
end, considering that the atoms composing each molecule conserve their identity, we 
expand the intermolecular potential energy as the sum of pairwise atomic energies. 
By expressing the intermolecular potential energy as a function of the intermolecular 
distance, we obtain a generalization of the usual multipolar expansion, containing 
additional terms scaling as inverse powers of the intermolecular distance. Those 
additional terms, which involve effective molecular multipole moments, are strongly 
anisotropic with respect to the molecular orientations. 

In the case of two vibrationally highly-excited Er 2 molecules, many additional 
terms bring a substantial contribution to the intermolecular potential energy. By 
estimating their characteristic lengths, we predict that those additional terms are 
also likely to influence the Er 2 -Er 2 collisions at ultralow energies. To conhrm that 
prediction, we can perform quantum-scattering calculations using the intermolecular 
potential-energy curves presented in this article. This would require however to know 
precisely the multi-channel wave function of the Feshbach states. Besides, since the 
intermolecular potential energy is a few mK (see Figj3]), which is the typical spacing 
between neighboring Feshbach levels, we can expect the intermolecular interaction to 
couple different Feshbach levels, and so to induce inelastic collisions. 
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Appendix A. Second-order correction and irreducible tensors 


In order to write Eq. (flTD as a sum of irreducible tensors, hrstly we expand the product 
of atomic multipole moments as 




tinrii 




+^A 




fcA=Ki-d'l <lA=-kA 


Q" 

~^{^i/ii)kAgA 


(A.l) 


where are the coupled atomic multipole moments (see Eq. flT^ ). and similarly 

for atoms j and j' of molecule B. Then we apply the transformation (see Ref. [38], Ch. 8, 
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Eq. (20), p. 260) 

/3-teifi tr 


c b a 
f e d > , 
j 9 t _ 


(A. 2 ) 


where [xiX 2 ---Xr] = ( 2 xi + l){ 2 x 2 + 1 ) x ... x ( 2 x„ + 1 ) and the number between curly 
brackets is a Wigner 9-j symbol, to a = f, b = £*/, c = if, d = i, e = ii, f = ij, g = k^, 
j = ks and t = k and to the corresponding components. We used that a 9-j symbol is 
unchanged after a permutation of two rows followed by a permutation of two columns. 

Secondly we work out the effective molecular multipole moments. We expand the 
products of Racah spherical harmonics in Eq. ffTOj) as (see Ref. [38], Ch. 5, Eq. (9), 
p. 144) 
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(A.3) 


where, recalling that 0 ^ = ©i = tt — ©2 and = tt -|- $ 2 , we used the property 

Cx 2 ^J. 2 i^ 2 , ^ 2 ) = (“1)"^^ Cx 2 ii 2 {^a, “hA)- After writing a similar equation for molecule B, 
we apply again Eq. (lA.2p to a = X'—iiif = Xii + \f, b = Xii, c = \f, d = X — iij = Aj + Aj, 
e = Xi, f = Xj, g = ka, j = Kb, t = IT, and the corresponding components including 
T = a. Then we apply the formula (see Ref. [38], Ch. 8, Eq. (26), p. 261) 


E faa fdS nby nc-i _ ftr naa 

^b/Sc-f^eefip^eegri^fipjlJ, [OeUL\ / ngi 




tr 



(A.4) 


to a = k, b = iij, c = ii'j>, d = 71 , e = Xi-\- Xj, f = A*/ -|- Ay, g = X, j = X', t = k, and the 
corresponding components including r = g. Then the invariance of the 9-j symbol: (i) 
by reflection about the anti-diagonal; (ii) by permutation of the resulting first two lines 
and the hrst two columns. Afterwards we deal with the Racah spherical harmonics for 
the intermolecular axis. Applying (see Ref. [3H], Ch. 5, Eq. ( 10 ), p. 144) 


Ec2wC„„(e, f = cioCc(«. 

aj^ 


to a = ka, b = Kb, c= k, a = qa, fd = Qb and j = g, we get to Eq. (ITHil . 


(A.5) 
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